Immunization Dynamics on a 2-layer Network Model 
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We introduce a 2-layer network model for the study of the immunization dynamics in epidemics. 
Spreading of an epidemic is modeled as an excitatory process in a small-world network (body layer) 
while immunization by prevention for the disease as a dynamic process in a scale-free network (head 
layer). It is shown that prevention indeed turns periodic rages of an epidemic into small fluctuation. 
The study also reveals that, in a certain situation, prevention actually plays an adverse role and 
helps the disease survive. We argue that the presence of two different characteristic time scales 
contributes to the immunization dynamics observed. 

PACS numbers: 87.19.Xx, 89.75,Hc, 87.23.Ge, 05.65.+b 



Recently the agent of severe acute respiratory syn- 
drome (SARS), a novel coronavirus, has spread over Asia. 
The epidemic apparently originated in the Guandong 
province of the People's Republic of China, and then 
spread rapidly throughout the world via air travel Q. 
For physicists the dynamics of epidemic spreading and 
immunization has been one of the fairly intriguing sub- 
jects. An amount of epidemiological work has been done 
on various networks; epidemic models of an infectious 
disease on small-world networks and computer virus 
infections on both original and structured scale-free net- 
works 0, Q , which seems to be due to the topology of 
Internet, etc. The effect of immunization is also studied, 
which shows that the targeted immunization strategy is 
more effective than the random immunization to lower 
the vulnerability in complex networks 0E1 In these pre- 
vious Refs. the effect of immunization on selected nodes 
is permanent and initially vaccinated nodes cannot infect 
or be infected by others in the future so that they func- 
tion as blocked passage or are merely removed from the 
viewpoint of the disease. Therefore this immunization 
by prevention is concerned with the static network struc- 
ture rather than some dynamical properties of models. 
For the more realistic prevention models let us consider 
the prevention as dynamical as the disease. 

Various information, which covers all the interactions 
among the elements of a system, spreads through the 
complex network. In a social network, one usually ob- 
tains many kinds of information through more than one 
channel, and which one is employed depends on the char- 
acteristics of the information, e.g. local interactions by 
physical contacts vs. global ones by mass media or In- 
ternet, public ones by way of official announcements vs. 
private ones at the club, etc. One can separate sev- 
eral kinds of channels from the whole network and make 
multi-layered network model, which means that the dif- 
ferent kinds of information propagate through the differ- 
ent layers and interact with each other when they meet 
at a node. The hybridization of networks was already 
attempted in Ref. jjj, indicating that the regular lat- 
tice superposed with the sparse random graph shows the 
clustering among the nodes and the small-world effect si- 



multaneously. Here we propose a 2-layer network model 
composed of a small-world network and a scale-free net- 
work, then analyze how the different kinds of informa- 
tion, which are the disease spreading and prevention in 
this study, compete with each other. 

Small-world and scale-free networks are the most fre- 
quently observed types, and which one is more proper 
again depends on what we are concerned about. For ex- 
ample, it can be a scale-free network in case of a venereal 
disease since it is found that the number of sexual part- 
ners shows a scale- free distribution || . Generally speak- 
ing, however, interactions between individuals arc limited 
by physical restriction such as spatial localization and the 
cost of adding social links 9] , henceforth most social net- 
works would be modeled by the single-scale network, e.g. 
Watts-Strogatz small-world network, or multi-scale net- 
work rather than a scale-free one. On the other hand, 
we expect a scale-free network to be suitable for spread- 
ing news or for governing the whole system owing to the 
existence of hubs and hierarchical structure. Since the 
prevention does not usually occur by physical contacts 
but most by means of communicating media or by some 
health policy, we separate a prevention strategy using a 
scale-free network, such as Internet, where an infected 
individual is supposed to report the status to his or her 
on-line nearest neighbors. Thus the existing prevention 
strategies should be modified into scale-free ones. We 
investigate a seemingly common situation where we con- 
trol the small-world excitatory system through scale-free 
inhibitory system. 

We prepare a ground where an epidemic and immu- 
nization by prevention directly compete with each other, 
so that we may get more information on the dynam- 
ics. For that purpose we introduce a 2-layer epidemic 
spreading and immunization model, or briefly '2-layer 
model'. Figure 1 displays the schematic structure of the 
model; the lower layer (body layer; BL) is a physical net- 
work over which the infectious disease spreads, while the 
upper one (head layer; HL) is an information network 
used for immunization by prevention. We notice that 
each channel has a different connectivity structure. We 
adopt here the well known susceptible- infected-refractory 
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FIG. 1: The epidemic spreading and immunization model 
with 2-layer structure. The infectious disease spreads over 
the Body Layer while the information for prevention does over 
the Head Layer. A susceptible node evolves to get infected 
or immune by its infected neighbors both on the HL and the 
BL. 
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FIG. 2: The infection and prevention cycle of a node. At 
once a susceptible node becomes infected it goes through the 
infection cycle deterministically. When it has /-state it can 
contaminate its neighbors. After the infection period it be- 
comes refractory and during the recovery period it becomes 
immune but not contagious. Prevention by its infected neigh- 
bors lasts for a prevention period. If nodes are immunized by 
a biological method tr and rp must be the same. 



(SIR) model [lOj for the characteristic states of a disease 
so that the health of each node is described by one of S, 
I, and R. An S-state node, when infected through con- 
tagion, starts an infection cycle of disease, i.e. it passes 
to the /-state, where it remains for a certain infection pe- 
riod Tj and then to a recovery period Tp before returning 
to the S-state. If an S-state node gets immune through 
the information channel, it stays for a prevention period 
Tp, then returns to the S-state, which is a prevention 
cycle (Fig. 2). An individual can get infected only dur- 
ing the S-state and only by contacting /-state neighbors, 
while it neither infects nor is infected by others during 
the Zi-state. In this sense /?-state nodes in both cycles do 
the same function as immune ones. Therefore it is nec- 
essary to divide the immunization into two classes; one 
is through infection and the other is by prevention. The 
existing prevention strategies correspond to the case of 
assigning the infinite Tp for initially selected nodes and 
with no further immunization. For a dynamical immu- 
nization we set both Tp and Tp to be finite and there 
is no initial selection of immune nodes. If the S-state 
nodes are immunized only by the immediate inoculation, 
Tp and Tp must be the same. Now we complete the 2- 
layer SIRS spreading and SRS prevention model. 

In the original small-world model, the rewiring proba- 
bility p characterizes the degree of disorder, ranging from 
a regular lattice to a completely random graph [Tl], [lj| , 
and is defined here as the number of randomly rewired 
links divided by the total number of links of the initial 
regular lattice. Hence at p = 0, N nodes constitutes a 
simple regular lattice each node of which has ko nearest 
neighbors, and becomes more random as p increases. We 
choose ko large enough to guarantee that the network 
shall be connected even after many random rewirings. 
Even at low p, the network shows small-world effect, i.e. 
the characteristic path length L(p), the averaged mini- 
mum path length between two randomly chosen nodes, 
increases logarithmically with the system size, while L(j>) 



does linearly in case of the purely regular lattice. On the 
other hand, a scale-free network is constructed by two 
rules: growth and preferential attachment 0,0]- Start- 
ing with a small number (too) of nodes, we add a new 
node with m links at every time step. Then the proba- 
bility that a new node is connected to node i is directly 
proportional to fej, the number of neighbors of node i. 
The resulting probability distribution of the node con- 
nectivity, P(k), follows a power law: P(k) ~ k~ 3 . The 
power law implies rather frequent emergence of highly 
connected nodes, often called hubs, which drastically re- 
duce path lengths. 

Let us define necessary parameters for the model study. 
A node with S-state may have kj, infected neighbors, 
each of which has probability q\> of contagion. At the 
same time, the act of prevention triggered by kh infected 
neighbors, each of which has the probability of im- 
munization by prevention, protects a node from being 
infected for a prevention period. With these provisions, 
a node is expected to evolve according to the following 
probabilities: 
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Starting with mo = to = 3, we construct HL as a 
Barabasi- Albert scale-free network with N — 10 4 and BL 
as a Watts-Strogatz small-world network with N = 10 4 , 
ko = 6, and p = 0.5. The infection period and the recov- 
ery period (= the prevention period) are set to be 4 and 
9 time steps, respectively, and initial infected fraction is 
to be 0.1. These values are selected following Ref. 0. 
Other initial conditions on various combinations of net- 
works have been also explored as well, and no significant 
changes are observed. 

We have performed a simulation with various values 
of probabilities and qi,. The density of the infected 
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nodes of the system, denoted by <Ji(t), is an important 
dynamical variable in this study. Figure 3 shows the evo- 
lution of <Ji(t) for various values of qh and <#,. The results 
are classified into three categories depending on infection 
probability qb- (1) For small values of qb, it is found that 
<Ji(t) vanishes in a finite time for all range of qh (Fig. 
3(a)). Here it is found that immunization by prevention 
lowers the peak of ai(t). (2) The results for an intermedi- 
ate range of qb are shown in Fig. 3(b). When prevention 
is absent for this range, Fig. 3(b) shows that cri(t) glob- 
ally oscillates in the same way as was reported in Ref . . 
When prevention is in effect, the amplitude of the global 
oscillation gets smaller as the prevention probability qh 
strengthens. (3) For large values of infection probability 
qb = 0.5, we observe a rather unexpected result. Figure 
3(c) displays that prevention plays an adverse role in this 
case and actually helps the disease survive, rather than 
eradicate it. 

For a discussion, it is useful to introduce two character- 
istic time scales; one is the characteristic time scale of the 
disease (Td ~ tj + tr) and the other is the time taken for 
the disease to sweep over the whole system (T s ), which 
is a function of qb, qh and we write its simplest form as 
Lb(^r^+T~Rqh), where Lb is the characteristic path length 
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FIG. 3: The numerical simulations where the BL is described 
by small-world network and the HL is described by scale-free 
network. A solid line represents qh — (no prevention), and 
a dotted line does qh = 1 (perfect prevention), (a) For small 
value of qb, (b) for an intermediate range of qb, where the 
dashed line is for the case of qh = 0.04, (c) for large value of 
Qb- 



of the BL. If Td is translated into the consumption of a 
species, T s would be the natural resources for the species. 
In case of T s ~ Td, let us assume that initially infected 
nodes are distributed locally, then the disease spreads 
over the whole system. As the almost infected nodes end 
their infection cycles and few susceptible nodes are left in 
some time, <Ji{t) shows a steep collapse (Fig. 4(a)). This 
is the way the system exhibits oscillatory behavior, and 
its period is expected to be slightly longer than Td- On 
the other hand, when T s is much smaller than Td, then a 
disease soon cannot find any more S"s left and becomes 
exterminated. Infection occurs only when /-state node 
meets 5-state node and is obstructed by i?-state node. 
Since qb is quite large, the disease spreads very rapidly 
over the whole system, reaching i?-state everywhere at 
some time, then we see that the disease cannot survive 
any more (Fig. 4(b)). Prevention, however, hinders this 
situation by always preserving S'-state nodes, and thus 
makes the disease a chronic one in the end. This does 
resolve the paradoxical result in Fig. 3(c). 

Figure 5 shows a phase diagram in the (qh, qb) plane, 
where darkened areas represent the regions of eradication 
of the disease, while white areas represent the regions of 
survival. In the lower dark region disease is eradicated 
because it takes too long time to prevail, i.e. T s 3> T&- 
On the contrary, it also does when it prevails too fast for 
individuals to recover themselves, as in the upper dark 
region where T s <C X^. The white regions corresponds to 
the situation when T s ~ Td- For qh = 0, one can find 
the range T s q < T s < T s i where the disease survives. 
The survival region in the plane of (qb, qh) in Fig. 5 fits 
well to this range of T s . For low values of qh, preven- 
tion becomes significant for the regions where qb is high, 
but becomes insignificant in the region where qb is low. 
Some analytic explanation would be helpful. Assuming 
that the system is continuous in time and each S-state 
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FIG. 4: Spatial characteristics showing the competition of two 
time scales. Underlying networks are 2-dimensional small- 
world networks with N = 50 x 50, k = 4, and p — 0.1. 
(a) The speed of spreading is relatively low, the oscillatory 
behavior establishes, tj = 5, tr = 10, and qb = 0.35. From 
the definitions of two time scales, Td ~ 15 and T s ~ 4.66. (b) 
The disease spreads over the system before initially infected 
nodes become refractory, tj = 7, tr — 14, and (ft = 0.5, thus 
T d ~ 21 and T s ~ 2.33. 
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node with k neighbors has the probability to have kai in- 
fected neighbors, we propose the following three ordinary 



d<J S (t) 
dt 

da jit) 
dt 

d<JR{t) 
dt 



where as, 07 and or represent the densities of sus- 
ceptible, infected and refractory nodes, respectively, pu 
and p' k are the normalized connectivity distributions of 
BL and HL networks, and / and g are monotonically 
increasing functions of an and 07 respectively, satisfy- 
ing f(0;Tn) = g(0;ri) = 0. tr and 77 mean relaxation 
times. Note that this set of equations already satisfies 
vs{t)+ai{t) + a R {t) = 1. 

From the stationary condition, one can find a thresh- 
old value of qb as a function of the average connectivity 
(k) = ^2 k kpk and the infection period 77, above which 
the disease survives. Assuming 5(07; 77) = ^j-, in the lim- 
its of <t 5 — ► 1 and a-/, 0ij — > 0, we obtain the following 
analytic results for the threshold value of 

9 Mfl = l-e- 1/(k)T '. (3) 
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FIG. 5: The phase diagram which shows the regions of erad- 
ication (dark) and survival (clear) activity of the system in 
the (qu, qb) plane. Each point is averaged over 100 realiza- 
tions. For low value of qt, the immunization by prevention is 
irrelevant. For low value of qu the effect of prevention appears 
clearly. 



differential equations in general: 



(2) 



The fact that qb t th is n °t a function of qu explains why 
immunization by prevention seems almost irrelevant in 
the low qb region. In Fig. 6, the analytic results are 
compared with the numerical calculations, ft becomes 
clear now that the shorter the infection period, the larger 
values of qb,th are needed for the disease to survive. 

Lastly we investigate the effect of variation of tr on the 
levels of small fluctuation 07 j in case of perfect preven- 
tion, qh = 1, in Fig. 3(b). For larger tr the levels of fluc- 
tuation is expected to lower, which coincides with the nu- 
merical simulations. Figure 7 shows the Boltzmann sig- 
moidal behavior of the relation: aij = T ^_ B + C, where 
A, B, and C are the positive constants, and a w 0.77 
from the fitting. For large tr the equation approximates 
to the power law with the exponent —a. We can explain 
the above result from Eq. (2). In case of qh = 1 three 
variables have the nonzero fixed values which are asj, 
aij, and <tr,/- From the stationary condition, assuming 
g(ar,Ti) = ^ again, a T j = {f(a RJ ;T R ) - as.^n is ob- 
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FIG. 6: Threshold values above which the disease survives. 
The curves show the equations for qt,th as a function of rj. 
The diamonds, circles and squares show the q bit h in the case 
of (k) = 4, 6, 8 respectively by numerical results. 
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FIG. 7: The levels of small fluctuation ajj with perfect pre- 
vention, qh = 1, in Fig. 3(b), as a function of the recovery 
period tr. The relation shows the Boltzmann sigmoidal be- 
havior. The dots represent the numerical results and the line 
does the sigmoidal fitting. 



tained. /(ctr./iTr) ~ since tr means the relaxation 
time. Even for the infinite recovery period the disease 
does not eradicate because C ^ 0. 

To summarize, we studied the immunization dynam- 
ics in epidemics by introducing a 2-layer network model, 
where an epidemic spreads on a small-world network of 
the system while the dynamical immunization by pre- 
vention takes place on a scale-free network of the same 
system. The prevention process indeed turns the global 
periodic rages of an epidemic into a small fluctuation. It 
is also found that prevention may help the disease sur- 
vive when certain conditions are met. We introduced two 
time scales, T s , Td, and explained the numerical results. 
Recalling that we did not include any genetic processes 
here, such as occurrence of immunity-resistant disease 
through mutation, all of these should be purely dynami- 
cal effects. 
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